Ferro-Nematic ground state of the dilute dipolar Fermi gas 
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It is shown that a homogeneous two-component Fermi gas with (long range) dipolar and short- 
range isotropic interactions has a ferro-nematic phase for suitable values of the dipolar and short- 
range coupling constants. The ferro-nematic phase is characterized by having a non-zero magneti- 
zation and long range orientational uniaxial order. The Fermi surface of spin up (down) component 
is elongated (compressed) along the direction of the magnetization. 
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Cold dipolar Fermi gases have attracted much atten- 
tion due to the novel anisotropic and long-range charac- 
ter of dipole-dipole interactions. Recent studies of many- 
body effects predicted an elongated Fermi surface (FS) in 
a one-component fully polarized Fermi gas with dipolar 
interactions along the polarization direction established 
by an external field [1, 2]. A biaxial state with a criti- 
cal value of the effective coupling constant = np? /cf 
was proposed in Ref. [3] (/i is the dipole moment of 
the fermion, n is the total density and ep is the Fermi 
energy of the free Fermi gas at the same density). It 
was found[3] that the system will exhibit violations of 
the Landau theory of the Fermi liquid both at quantum 
criticality and in the biaxial phase. More generally, un- 
derstanding anisotropic non-Fermi liquid phases of cold 
atomic systems may shed light into the quantum liquid 
crystal phases in strongly correlated systems and high T c 
superconductors [4-6] . 

The question we want to address here is whether a 
cold spin-Y2 Fermi gas with long range dipolar interac- 
tions can become spontaneously polarized and what is 
the nature of the broken symmetry state. The theory 
we present here is a generalization of the theory of the 
Stoner (ferromagnetic) transition in metals [7] to take into 
account the effects of the long range and anisotropic dipo- 
lar interaction [8]. As we will see below, much as in the 
theory of Stoner ferromagnetism, the polarized state can 
occur only for sufficiently large values of the magnetic 
dipole moment and/or of the spin-flip scattering rate. 
However, unlike what happens in Stoner ferromagnetism, 
as a result of the structure of the dipolar interactions, 
the resulting polarized state is also spatially anisotropic, 
a ferro-nematic state. 

The classical version of this problem has been con- 
sidered in mixtures of ferromagnetic particles with ne- 
matic liquid crystals[10], and in dipolar colloidal fluids 
and ferrofluids[ll]. Classical dipolar fluids have com- 
plex phase diagrams, typically featuring inhomogeneous 
phases with complex spatial structures. Much less is 
known about their quantum counterparts. In the case 
of simple quantum fluids, such as 3 He, the dipolar inter- 
action plays a small role compared to the short-range ex- 
change interaction[12]. In the context of ultracold gases, 



a number of atomic and molecular systems with strong 
dipolar interactions, such as Dy, have been the focus of 
recent experiments (see Ref. [13]). 

Consider a restricted Hilbert space of two hyperfine 
states, called 1 ("spin up") and 2 ("spin down"), of a 
point-like magnetic atom of mass m and magnetic mo- 
ment, M, with components Mi = fia 1 (i — 1, 2, 3) and a 
are the usual spin- 1 /^ Pauli matrices (the factor of H/2 is 
absorbed in the definition of fi) . The Hamiltonian is 

H = J d 3 x^(x) (-^") VA*(x) + ~ J d 3 xd 3 x' 

, X% Q ';/3 / 3<V>/3'( X ')Vv(x) (1) 

where the fields Va( x ) destroy fermions on spin state 
with z-component a — 1,2 at position x. We consider 
the model interaction is of the form 

a 2 . 

[/ (x,x) aaW = —a^idij - 3Tifj)a$p, 

+ 9 <W<W<5(r) (2) 

where r = (x — x')/|x — x'| and f is a unit vector in 
the direction of r. The last (ultra-local) term represents 
the short-range isotropic (contact) interactions. It only 
affects the spin-triplet channel and we denote by g the 
associated coupling constant. The Fourier transform of 
the bare two body interaction is (47r/z 2 /3)cr^ ct / (3qiqj — 
5ij)o~Lp, + gSaa'dpp' . The Hamiltonian is invariant under 
simultaneous SU (2) transformations in spin space and 
50(3) rotations in real space. This mixing of orbital and 
spin degrees of freedom is, in essence, what relates the 
distortions of the shape of FS and the spin polarization. 

The ferro-nematic state breaks simultaneously the ro- 
tational invariance in spin space and in real space of the 
Hamiltonian, and its order parameters reflect this pattern 
of symmetry breaking. The order parameters are: a) the 
local magnetization vector M a (a = x, y, z) that measures 
the spin polarization, b) the nematic order parameter, a 
3x3 symmetric traceless matrix, Qij (i,j = x, y, z) that 
measures the breaking of rotational invariance in space, 
and c) the generalized "nematic-spin-nematic" order pa- 
rameter Qij, a tensor symmetric and traceless on the 
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FIG. 1: (color online) Schematic phase diagram as a function 
of A s = gn/tF and Ad = nfi 2 /ef- It has a paramagnetic phase 
(M — 0), and a ferro-nematic phase with partial (0 < M < 1) 
and full polarization (M = 1). Dashed (full) curve: first 
(second) order phase transitions; A: tricritical point. 



spatial (i,j = x,y,z) and spin (a, b = x,y,z) compo- 
nents, that measures the breaking of both symmetries 
[14]. General order parameters of the latter (nematic- 
spin-nematic) type were considered by Wu et. al. [16] 
who gave a detailed description in 2D systems and par- 
tially in 3D systems. The ferro-nematic state has an 
unbroken uniaxial symmetry in real space. It is a 3D 
generalization of the quadrupolar a phase of Ref. [16]. 

An intuitive way to describe these phases (and their 
order parameters) in a Fermi system is in terms of spon- 
taneously deformed Fermi surfaces [3, 5, 16]. A ferromag- 
netic state is isotropic in real space and has a spherical 
FS of unequal size for both spin polarizations. The ne- 
matic phase is isotropic in spin space, and its FS has a 
uniaxial distortion in real space, with the up and down 
spin FS being identical in shape and size. For small val- 
ues of the order parameter, the distorted FS are ellipsoids 
with an eccentricity determined by the magnitude of the 
order parameter Qij . Due to the mixing of orbital and 
spin degrees of freedom by the dipolar interaction, fer- 
romagnetism causes the FS to distort, thus driving the 
system into an uniaxial nematic state. Since the state 
is ferromagnetic, the up (down) FS is a prolate (oblate) 
revolution ellipsoid. Both FS are collinear, and have un- 
equal distortion and volume. Hence, the nematic-spin- 
nematic order parameter has a finite value. We find 
that all three order orders are present even for arbitrary 
small values of the dipolar coupling, where the ferromag- 
netic is the strongest, the nematic intermediate, and the 
nematic-spin- nematic the weakest. 

In this paper we use a Hartree-Fock (HF) variational 
wave function to determine the phase diagram as a func- 
tion of the density n, the (dimensionless) dipolar coupling 
constant A^ and of the (dimensionless) local exchange 
coupling A s = gn/tp. We consider an infinite system 
and significant finite size effects, such as the trap poten- 
tial and the associated inhomogencity of the gas, are not 
considered but can be included using the Thomas-Fermi 
approximation. We find two phases: a) an isotropic un- 



polarized state and b) a ferro-nematic phase. As in the 
conventional theory of the Stoner transition, we find that 
the values of A s and A^ on the phase boundary are of or- 
der unity. In this regime, a HF wave function can only 
yield qualitative results, such as the broad structure of 
the phase diagram, but its is not not expected to be quan- 
titatively accurate. Even within mean field theory, HF 
may miss important physics; a recent study [17] found 
that to second order in g the usual continuous Stoner 
phase transition can turn first order. 

We take a variational HF wave function of the form of 
a Slater determinant describing a state in which the spin 
up and down Fermi surfaces are spontaneously deformed 
away from their non-interacting spherical shape. We will 
not consider other interesting states with more complex 
order, such as biaxial[3] and its generalizations. Since we 
are interested in magnetism we allow for the volume of 
the up and down Fermi surfaces to change as well. This 
results in a ("Thomas-Fermi" like) distribution function 
of fermions in momentum space with 4 variational pa- 
rameters, kpi, k>F2 oi\ and a 2 . We keep the total particle 
density n = n\ + n 2 fixed, [1, 2] 
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where Q(x) is the step function and a = 1,2. If a a 
I both FS's are spheres. Eq.(3) has the property that 



the total density does 
with 



not depend on the FS distortion parameters a 
a a > 1 (oblate), a a < 1 (prolate). 

Computing the energy in HF we obtain an expres- 
sion for the ground state energy density in terms of 
the distribution functions of spin up and down particles 
?Vk = (c^. k c CTk ), and J2k n <?.M = N v, 
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where #k-k' labels the direction of k — k' with respect 
to the z-axis. £j ni is a function of the magnetization 
AI = (ni — ri2)/n, the density n, and the ratio of the 
undistorted Fermi surface volumes is r 3 = k F1/ 
(1 + M)/(l - M). The total energy density E is 

dh 2 
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where C x = (3/10)(6tt 2 ) 2 / 3 . The integrals 1(a) and 
I(ct\, a>2, M) are given by 

T , % „ 6 3arccosa 3 / 2 , . 

1(a) = -2 3 — r- T -^ i 2 -r m 6 

a 6 — 1 (a L — a z ) A ' z 
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The total energy £ becomes 

V + M)«/»(2a 1 + ^) + (l- M )«/»(2« a + ^) 

(1 + M) 2 /(ax) - 2(1 - M 2 )I( ai ,a 2 ,M) + (1 - M) 2 I(a 2 ) 
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where £o = (2Ci?i 2 n 5 / 3 )/(2 5 / 3 m). Numerically minimiz- 
ing the energy we obtain the phase diagram in the cou- 
plings X s and Arf of Fig. 1. 

The phase diagram shown in Fig. 1 exhibits a ferro- 
nematic phase and a paramagnetic phase, separated by 
a phase boundary consisting of a line of 1st order tran- 
sitions that meets a line of continuous transitions at a 
tricritical point (labeled by A.) As expected the ferro- 
nematic state becomes more accessible as the s-wave cou- 
pling increases. This phase is fully polarized (M = 1) 
for most of the phase diagram, except for a small re- 
gion where the polarization is partial, < M < 1. In 
the ferro-nematic phase with partial polarization the up 
and down Fermi surfaces are unequally distorted, while 
in the fully polarized regime only one distorted up FS 
exists. Partial magnetization in the conventional Stoncr 
transition occurs for = and 4/3 < A s < 2 2 / 3 , where 
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FIG. 2: (color online) Magnetization M and Fermi surface 
distortion parameters ai and oli vs the dimensionless dipolar 
coupling X d for an s-wave coupling X B = 1.34. 




-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

&d ^d fcy 

(a)A 3 = 1.34, \ d = 0.04 (b)A s = 1.34, \ d = 0.08 




-2-1 1 2 -3-2-10123 



(c)A s = 1.37, A d = 0.10 (d)A s = 1.37, A d = 0.14 

FIG. 3: (color online) a-d: Spin up and down FS, labeled by 
ni and 712, for several values of A s and X d and fixed particle 
density n. a d — m/j, 2 /h 2 is the dipole-dipole interaction length 
scale. Dashed lines: Free fermion FS. 



the up and down FS become distorted even for arbitrar- 
ily small values of the dipolar coupling (see Fig. 2.) 
From the structure of the free energy for small values 
of Qi = cti — 1, Q2 = ct% — 1 and M (valid in the vicinity 
of the continuous transition) we see that the dipolar in- 
teraction leads to a leading term of the form M(Qi — Q2), 
which is invariant under 1 <-> 2. Since the ferromagnetic 
state is already favored by the contact term, this term 
also favors Q\ < and Q 2 > 0. 

In the ferro-nematic phase the SU(2) spin symmetry 
of the Hamiltonian is broken down to the residual U(l) = 
50(2) invariance of this uniaxial state. The equilibrium 
FS of the up and down spin components are shown in Fig. 
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FIG. 4: (color online) Pressure (broken), bulk compressibility 
(full) and chemical potential (dashed) (normalized to their 
values at the non-interacting Fermi gas) vs the dimension- 
less dipolar coupling A^, for X a — 0.8. At A^ ~ 0.3 there is 
a 1st order paramagnetic-ferro-nematic transition. The gas 
becomes unstable at Ad ~ 0.52. 



3 for several values of of the coupling constants. Both 
FS's are invariant under 5*0(2) rigid rotations about M. 

The total energy, the equilibrium values of the FS dis- 
tortions and the magnetization are functions of the par- 
ticle density n = N/V. The pressure P = —(dE/dV)N, 
chemical potential fi — (dE/dN)v and the compress- 
ibility K~ l = n(dP/dn) can be computed straightfor- 
wardly, and decrease monotonically as the dipolar cou- 
pling increases, see Fig. 4. The compressibility vanishes 
at A^ ~ 0.52 [2, 18] where the Fermi gas becomes (for- 
mally) unstable to collapse[19]. 

In current experiments on cold atoms it is possible to 
prepare a two-component Fermi gas out of the hyper- 
fine manifold of the atom. One can imagine, a two- 
component dipolar Fermi gas made out of, say, the 
hyperfme manifold of strongly magnetic atom such as 
Dy. For 163 Dy with a density of 10 1 cm" 3 , we estimate 
Xd — 0.01 [20, 21]. The recent experimental observation 
of itinerant (Stoner) ferromagnetism in ultra-cold gases 
of 6 Li atoms [23] opens the possibility to detect the ferro- 
nematic state in the laboratory possibly by tuning the 
s-wave scattering amplitude. Depending on the strength 
of the dipole moment and of the s-wave coupling, the FS 
for the two components may differ considerably from the 
spherical shape, and free expansion experiments may be 
able to provide signatures of this state. Ferromagnetism 
in cold dilute systems interacting only by s-wave contact 
pseudo-potential has been difficult to observe as the spin 
states of fermions are conserved separately on the time 
scales of the experiments [17]. However, a ferro-nematic 
state can be set up experimentally in systems with a pop- 
ulation imbalance of hyperfme states, which may also 
exhibit other anisotropic and inhomogeneous phases. In 
actual experiments the trap potential is anisotropic and 
weakly inhomogeneous. We estimate that such effects do 
not affect our main results provided the trap potential 
aspect ratios depart from 1 by perhaps up to about 30%. 



Nevertheless, the anisotropy acts as a weak symmetry 
breaking field, that orients the ferro-nematic order. 

In this work we have shown the existence a new phase 
of matter, the ferro-nematic Fermi fluid, a ground state of 
a dipolar Fermi gas with short range interactions with a 
spontaneous magnetization and long range orientational 
order. In this state, the up and down FS manifolds have 
unequal shapes and volumes. Since rotational invariance 
in real space and in spin space is simultaneously sponta- 
neously broken in this state, it supports a rich spectrum 
of Goldstone excitations. As a result the fluid is an opti- 
cally anisotropic medium whose properties that may be 
detected by light scattering experiments. 
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